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A scattering of short, weak, nonresonant laser pulses on a 
Bose-Einstein condensate is proposed as a tool for studying its 
statistical properties. We show in particular that the variance 
of the number of scattered photons may distinguish between 
the Poisson and microcanonical statistics. 
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I. INTRODUCTION 

The mean field approach to describe properties of 
Bose-Einstein condestates (BEC) in weakly interacting 
atomic gases Jl|-|| has proved to be immensly succesful 
B. However, quantum statistical properties of the con- 
densate, and in particular higher order correlations in 
BEC have been a subject of considerable interest in the 
recent years. In the series of papers the fiuctations of the 
number of condensed atoms of the ideal Bose gas were 
calculated in micro- and canonical ensembles The 
corrections due to the weak interatimic interactions were 
considered by several authors, but different models give 
different results, so that the impact of collisions on the 
statistics remains unclear [Q. The possible reason for 
these diffuculties is the ambiguity in defining condensed 
fraction. Similar ambiguities lead to difficulties in deter- 
minations of the modification of the critical temperature 
for weakly interacting gases || . 

While the mean number of condensed atoms Nq as a 
function of temperature has been measured in several lab- 
oratories J9| , from the experimental point of view very lit- 
tle is known about fluctuations of BEC. The main source 
of information concerning higher order correlations comes 
from the studies of the depletion of BEC due to inelas- 
tic two-body and three-body collisions JlO| . In this way 
the mean value of Nq and Nq, have been estimated [ plj . 
The experimental result have unequivocally ruled out the 
thermal fiuctations in the condensate. Precision of those 
measurements is, however, unable to distinguish between 
sub- and super-Poissonian, but normal fluctuations. 

In fact different models of the interacting atoms con- 
densate coexist in literature. In particular, the classic 
Bogolyubov approximation |l^] tends to favour the Pois- 
son statistics of the condensate. It would be desirable to 
discriminate between the models with the help of a clean 
experiment. 

This paper describes such a proposal. We propose to 
use a scattering of a short weak non-resonant laser pulse 



as a means of probing the BEC statistics. The problem of 
light scattering on the condensate has been thoroughly 
studied in the recent years Most of these papers, 

however, concerned the signatures of BEC transition in 
scattering, rather statistics of the condensate. The refer- 
ences most closely related to the present paper are papers 
by You et al. |l~4|,[l5| on scattering of short laser pulse 
on bosonic and fermionic clouds. These papers use very 
similar approximations to solve time-dependence of the 
system, but both use grand canonical ensemble, and do 
not address the problem of the condensate fluctuations. 

The paper is organized as follows. In section II we 
present the theoretical model, and solve it for the pho- 
ton operators of the scattered field. In Section III we 
analyze the mean number of scattered photons, whereas 
in Section IV its variance in function of the assumed sta- 
tistical properties of the condensate, and physical pa- 
rameters such as temperature and scattering angle. We 
compare the result for the thermally fluctuating gas, for 
the coherent state, for the ideal Bose gas described by 
the microcanonical ensemble and for the nonfluctuat- 
ing gas. While the mean number of scattered photons 
can only discriminate between the thermal and the con- 
sidered non-thermal condensate states, the variance can 
further distinguish the coherent state from the last two 
models. In Section V we present some conclusions. 



II. THE MODEL 

Our method of exploring the statistical properties of 
the condensate is based on the scattering of the series 
of short light pulses. The fluctuations in a given con- 
densate occur at the time scale given by the interatomic 
collisions, thus to probe the fluctuations the time de- 
lay between consecutive pulses should be of the order of 
miliseconds. The distribution of the number of photons 
scattered into a given solid angle should be measured. 
Out of the distribution of the number of scattered pho- 
tons we may compute the mean and its variance. We 
assume that the pulse of light is weak and far detuned in 
order to avoid heating the condensate during the interac- 
tion. The pulse of light should be also sufficiently short 
in time. It should satisfy two conditions: 1. uj rnax t <C 1 
where huj max is the energy of the highest occupied level 
in the trap, since we are going to ignore factors e luJnt , 2. 
t <C 7 _1 or the lifetime of the transition, since we are 
going to ignore the spontaneous emission. Note, that the 
conditions for the proposed experiment are complemen- 



1 



tary to the well known method of phase contrast used 
for nondestructive imaging of the condensate 16 1. In 
fact we are interested in the regime in which the (small) 
condensate does not lead to significant refraction, while 
the scattered light carries then more direct information 
about fluctuations. 

The total Hamiltonian consists of the following parts 



H — H a + H a i + H a f + H 



(1) 



The term H a is the atomic Hamiltonian. In the second 
quantization formalism it reads 



H a = Ku a n9n9fi 



-wo)ete,R, (2) 



where Uu)Z and hu^ are the energies of the CM mo- 
tion of atoms, in the ground and the excited electronic 
state respectively. The oj is the atomic resonance fre- 
quency. The operators g^ (e^) arc annihilation oper- 
ators of atoms in ground (excited) electronic state re- 
spectively. They fulfill bosonic commutation relations 
[9H,9a} = <5n,m, [e«, et] = <5s,m- We use vector indices n 
and m because we consider a three dimensional trap. 

The term Hf in the Hamiltonian is the energy of the 
electromagnetic (e-m) field: 



x 



H f = / d 3 k hujrat a r , 



k k \ kX> 



(3) 



where (^k.) are annihilation and creation operators 

of photons with wave-vector k (and frequency lo^ = kc) 
and polarization A respectively. 

The interaction of the atoms with the laser light is 
described by the H a i term: 

Hai = ^^{n,g\e- i ^\rn,e)e i ^ t gl t e rh + h.c. (4) 



where ujl is the laser frequency, f2 = 2Ed/h is the 
Rabi frequency with the transition dipole moment d and 
amplitude of the electric field of the laser E. In the 
above equation Franck-Condon factors (n, g\e~ lkLr \rh, e) 
appear. They are proportional to the amplitude of the 
transition between two states of the trap induced by ab- 
sorbtion or emission of a photon. The trapping potential 
felt by excited and ground-state atoms may in general 
be different, and therefore their eigenstates may also be 
different. The labels e and g distinguish these two sets of 
states. However, in later derivations we will assume the 
same potential for all atoms, and we will use notation |n) 
for the single atom trap states. 

The atom-field interaction part of the Hamiltonian has 
the form: 



H. 



k £(fj, g\e~ tkLp \m, e)gte m + h.c. (5) 



The Heisenberg equations of motion for the atomic op- 
erators gfi, can be solved in the following approxima- 
tion. The atoms are driven only by laser field dominat- 
ing over vacuum modes. In our approach we neglect the 
back action of atoms on the driving light and the influ- 
ence of the vacuum modes. This last approximation is 
just a neglect of the spontaneous emission. It is justi- 
fied if the duration of the pulse is short compared to the 
spontaneous emission life time. The laser light, as one 
may notice from H a i, is treated clasically. The equation 
of motion for atomic operators in the interaction picture 
with respect to H a are given by 

9n(t) = -i^e^^^ty^^iftle^A), (6) 
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where A = u>l — loq is the detuning of the laser. Now 
the assumption of the short pulse can be invoked to 
approximate the factors: e ! ^™~""- ) * « 1. Then, the 
system of equations @-(0) is easy to solve analyti- 
cally. To this end we introduce the operator Eft(t) = 



Em 1 



-ik L r 



m)e^(t). The system of infitely many cou- 



pled equations then reduces to a set with only pairwise 
coupling. In the next step, we insert the found solution 
into the evolution equation for operator ag, . In our ap- 
proach the atom charges and currents distribution is a 
source for the vacuum modes of the e-m field, while the 
emmitted photons have no influence on the evolution of 
atoms. In other words we neglect the interaction of scat- 
tered photons with other atoms. The equation of motion 
is solved by simple integration in time, which leads to 



Hx + \l 2fe (27T) 3 



x e 



i(oJg-wz,)t/2 



/i(*)£4< 



i(k L -k)r 



m)gjj 



Mt) 



YA(rh\e^ k ^\n)g n 



(8) 



where ft' = yj A 2 + fl 2 , and fiit) are functions of 

time given in Appendix A. All the operators on the RHS 
of Eq. (||) without specified depedence on time are taken 
at the time t = 0. We will use this convention in the next 
sections. The dots in the brackets denote the remaining 
terms proportional to gLe^ and eLe^i. They do not 
give any contribution to the mean number and variance 
of photons. 
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III. MEAN NUMBER OF SCATTERED 
PHOTONS 

The basic statistical quantity we can construct from 
the solution (0) is the mean number of photons. We 
would like to know if it depends on the statistical prop- 
erties of the condensate. In particular we are interested 
in its relation to the condensate fluctuations 6 2 N . In 
an experiment the angular distribution of scattered pho- 
tons may be measured by scattering the series of short 
pulses, then calculating the mean from a number of de- 
tected photons in every pulse at the given angle. The 
calculated value will be approximately the mean number 
of photons of given mode, summed over polarization and 
integrated over the frequencies: 



dN 



ph 



dn 



(0,<M)=£ Jdk k 2 {d\ x {t)d ix {t)). 



(9) 



The k under the mean value on the RHS of Eq (||) has a 
direction given by (9, </>). As we show later the spectrum 
has two spectral components, and we do the integration 
in Eq. (^|) for every component separately. Hence, the 
integration leads to the total intensity of a given peak in 
the spectrum. The mean value under the integral should 
be calculated by tracing with the particular density ma- 
trix. Since we work in the Heisenberg picture, this den- 
sity matrix should be the initial one for the whole system: 
atoms and e-m field. At t = the density matrix is just 
the product of the density matrices of every subsystem: 
p(Q) = p g {0) <g> j0 e (O) ® /Q/(0), where letters g, e, / cor- 
respond to ground state atoms, excited atoms, and E-M 
field respectively. We can assume that for t = all atoms 
were not excited, hence p e (0) — |0, 0, . . .)(0, 0, . . . | 

Our approach pertains to a so called weak condensa- 
tion regime, that is the situation in which the number of 
atoms in the trap is not too large, and the fluctuations are 
therefore relatively more important. In this situation the 
role of interatomic interactions is to assure the migration 
of the system in the phase space, so that the density ma- 
trix can be well described by a density matrix of an ideal 
gas in microcanonical, or alternatively canonical ensem- 
ble. Experimental conditions favour the microcanonical 
description of the condensate (energy and the number of 
particles conserved), therefore the density matrix p g (0) 
is determined by the microcanonical ensemble Q . At the 
initial moment the only occupied modes of the e-m field 
are these corresponding to the laser light. But the occu- 
pation of these particular modes does not affect the re- 
sults of the measurements that are done for the nonzero 
angles. Therefore we put p/(0) = |0, 0, . . .) (0, 0, . . . | at 
the initial moment. 

Substituting the annihilation and creation operators in 
Eq. (|J) by explicit formula (|J), after some calculations 
(see details in Appendix A) we obtain 



dN 



pit 



d 2 n 2 



dn dt 



32ir 2 e frc 3 A 2 



{\-{n k n d ) 2 )^f{q), (10) 



and 

(ii) 

where n k is the unit wector in the direction (#, <j>) and 
fid is the unit wector in the direction of the dipole mo- 
ment d. The frequency uj is the frequency of the given 
component of the spectrum. Following the calculation in 
Appendix A, the spectrum consists of two peaks with fre- 
quencies: luq (inelastic component) and lul (elastic com- 
ponent). The vector q is the vector of the momentum 
transfer: q = k — k^. The intensities of both elastic and 
inelastic components are almost the same. Note that in 
the presence of spontaneous emission losses the inelastic 
peak would be suppresed for longer pulses. 

For the numerical purposes our expressions ( pd| ) were 
modified. After some calculations the mean value in- 
volving the operators g may be replaced by the statisti- 
cal quantities for the microcanonical ensemble (see Eq. 
(B4) in Appendix B). The statistical moments can be 
expressed by the microcanonical partition function, and 
then calculated by means of the recurrence algorithms. 
In order to check if the influence of the fluctuations may 
be noticed in the measurments of scattered photons, we 
compare the microcanonical results to the scattering on 
the nonfluctuating condensate. Such state may be the- 
oretically realized in the following way: the mean oc- 
cupation number for each level is the same as for the 
microcanonical condensate, while the higher statistical 
moments are decorrelated 



(NiNj) = (Ni)(N 



3 1 1 



(12) 



where indices i and j label the states in the trap[]. From 
this property it directly follows that the condensate does 
not fluctuate: SNq = in this particular state. 

There are other possible statistical properties of the 
condensate. The time-honored Bogolyubov approxima- 
tion, taken at its face- value, assumes that the condensate 
is in a coherent state 113] . We therefore consider also the 
Poissonian distribution: 



(N 2 ) = (N) 2 + (N), 
(N l N 3 ) = (N l )(N J ), 



(13) 
(14) 



The assumption that excited states are also coherent may 
be questionable, but we do the calculations for low tem- 
peratures, where the statistical properties of the conden- 
sate are dominant. For completness of our review of the 



1 We omit the vector indicies because the mean ooccupation 
number depends only on the energy of each state and not on 
its spatial properties. Therefore the energy is sufficient to 
enumerate the states. 
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different states we included also the results for the con- 
densate described by thermal state. Of course this statis- 
tics, implicit in the grand canonical ensemble has already 
been ruled-out by the experiment jnj. The statistical 
properties of this last state are imposed by 



(Nf) 
(NiNj) 



2(N,) 2 + 



(Ni), 



(15) 
(16) 



The comparison of the scattering of photons on these dif- 
ferent condensates consisting of 1000 atoms is presented 
in Fig.l. We plot here the elastic part of the scattering. 
The conclusions for the inelastic component are very sim- 
ilar. The quantity proportional to the mean number of 
photons is plotted as a function of the angle 9. This mean 
photon falls-off exponentially for large angles due to the 
presence of the, so called, Debye factor (see Appendix A). 
The results were computed for realistic parameters. The 
oscilator length £ was calculated for the sodium atoms in 
the spherically symmetric trap with u> = 400s~ 1 . The 
resonance frequency of the transition loq corresponds to 
the wavelength Ao = 800nm, and detuning A = 1GHz. 
Analysis of Fig. 1 leads to the conclusion, that only scat- 
tering on the "thermal" condensate may be clearly dis- 
tuingished. There are, of course, differences between the 
scattering on the remaining three models, but they are 
small, and become even smaller if the number of atoms 
increases. Our main conclusion in this section is that 
the mean values do not allow to detect and measure the 
condensate fluctuations. 
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FIG. 1. Angular distribution of the mean number of the 
scattered photons measured with respect to the direction of 
the incident laser pulse in the plane perpendicular to the po- 
larization axis of the incident light. The elastic component of 
the spectrum is shown. Of four different density matrices de- 
scribed in the text only the result for the thermally fluctuating 
source may be distinguished. The mean number of photons 
is scaled by the corresponding results for a single atom. 



IV. THE VARIANCE 

As we have seen above, the mean number of scattered 
photons is fairly insensitive to the statistics of the con- 



densate. We have to look at higher moments of the pho- 
ton statistics. Its variance is a suitable quantity. It may 
be experimetally measured. It is sufficient only to scatter 
a series of pulses and after that to calculate the variance 
of the number of photons detected during every pulse at 
given angle. Repeating the series of measurements for 
different angles leads to the angular dependence of the 
variance. The formula which corresponds to this mea- 
sured quantity is the following 



X J 

\ 2 



(17) 



The mean values under the integrals should be calculated 
as in the previous case by tracing with density matrix 
of the initial state. The integration over the frequency 
is performed for each constituent of the spectrum sepa- 
rately, because we are interested in checking the variance 
of the photons both in elastic and inelastic scattering. 
Substituting the solution (^|) for the creation and annihi- 
lation operators in (|l7]), after a tedious but straightfor- 
ward calculations we obtain the following result 



\ dCldt J 



327r 2 e ?ic 3 A 2 



(l - (n k n d ) 2 ) y/g(tf), 
(18) 



and 



n,n\rh,fh r l,l f ,j ,j f 

x ( (9a9a> 9^9™* 9f9p gtgj, ) 

- (grJn'gl-Jrn')(glgpgj9f)) , (19) 

where we denote Frank-Condon coefficients by /3- -, (q) = 
(n\e^ p \n'). The Eq. @ is very similar to Eq.@ for 
the mean value, and we use the same notation. The 
frequency u> is, as previously, the frequency of the con- 
sidered component of the spectrum. In order to perform 
the numerical calculations the mean values of the atomic 
operators in formula (^) are expressed by means of the 
different statistical moments of the occupation number of 
the trap states. Unfortunately even when we couple the 
annihilation operators g with the creation operators g' 
into pairs, we still have to sum over four indices. These 
indices are additionally the vector ones. Furthermore the 
number of possible combinations of these pairs is much 
larger than in the case of the four operators in Eq.(|Tl|) 
for the mean number. We have to use some approxima- 
tion to be able to perform the summation for reasonable 
size condensates (~ 1000 atoms). We may utilise the fact 
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that for BEC the number of atoms in the ground state 
is of the order of the total number of atoms, so we ex- 
pand our formula ([l9|) in the powers of g^go- As it turns 
out, for 1000 atoms a sufficient approximation requires 
inclusion of fourth-, third- and second order terms in this 
expansion. In this case summation is performed over two 
different excited states, what is possible to realize. 

In Fig. 2 we present the variance of the scattered pho- 
tons calculated on the condensates with different statis- 
tical properties, described in the previous section. The 
numerical computations were realized for the same phys- 
ical parameters. Again we only present the results for 
the elastic component. For vertical axis we also choose 
the same units as in the previous picture. Now we are 
able to distinguish the scattering on the coherent state, 
what was not possible on the basis of mean numbers. 
The larger difference occures for the elastic scattering 
for small angles. The difference between scattering on 
the microcanonical and uncorrelated condensated is still 
very small. However, for small angles these curves are 
different. In this regime our approximation is not valid 
for the microcanonical condensate. The expansion of the 
variance in the powers of q, shows that the result be- 
haves as 9 2 and the curve is indistinguishable from the 
curve for uncorrelated condensate. As we expected, the 
scattering on "thermal" condensate leads to completely 
different results. In fact while microcanonical and un- 
correlated states are even harder to distinguish for larger 
number of atoms, the distance to the results for the co- 
herent states grows linearly with N. 
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FIG. 2. Angular distribution of the variance of the number 
of scattered photons is shown. As in Figjl| the elastic compo- 
nent of the spectrum is shown. Note that now not only the 
thermal but also the coherent states results are clearly distin- 
guish. It is still impossible, however, to distinguish between 
the microcanonical and the decorrelated results. 



V. CONCLUSIONS 

The statistical properties of the Bose-Einstein conden- 
sate are of great interest. Its measurments would advance 
the theory of the coherence of existing atom lasers. 



In the present paper we propose to use the scattering 
of short laser pulses for the study of fluctuations. One 
needs higher moments of the distribution of the number 
of the scattered photons to have a sufficiently subtle tool. 

The main open question is the influence of weak in- 
teractions on the statistics. While the theory remains 
muddled, we know that the main impact of interactions 
is the emergence of two length-scales: that of the size 
of the condensate, broadened by the repulsive forces and 
that of the size of thermal cloud, expelled from the center 
of the trap by the condensate. It means that we would 
have two different Debye factors damping the scattering 
cross-section at large angles. This way by changing the 
angle one could, perhaps, gain access to the fluctuations 
of the condensate and fluctuations of the thermal cloud 
separately. They are lumped together for the ideal Bosc 
gas studied in the present paper. 
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APPENDIX A: THE PHOTON SPECTRUM 

The function of time which appears in the solution (j^) 
are given by 



/i(*) = 



A 



. UJ k - UJ L 
■ sm( 1) 



h{t) 



UJ k ~ 0J L 

O'-A _m± 
e 2 

. .UJk- UJL j\ 

sm( 1) 

UJk - UJL 2 

_ in't 

~A 6 2 

4w 2 fc 



— e 2 sm(wifct) 

4wi fc 



sin(w 2 fc£), 



sin(w2fc£), 



±jj 



lk 



(Al) 

iil't . . . 

e 2 sm(ujikt) 

(A2) 



where we use the following notation: il' = VA 2 + Vl 2 
uj\k = (uJk - UJL + fi')/ 2 ) W 2A: = (uJk - uJl - fi')/2. Our 
laser pulse is weak and far detuned and we may simplify 
the equations ( A1)- (|A2[ ). For these assumptions: A 3> O 
and f2' s» A for positive detuning. The functions of time 
fi and fi take the form 
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hit) = — sin — t sin — t , 



V 2 



V 2 



/a(*) = — sin (-^tj - — sin f 



2ujfco V 2 
' OJk + lu - 2lu l 



2(uj k +lu - 2io L ) 



(A3) 



(A4) 



where tOkL = — and uj k Q — uj k — luq. Calculation 
of the mean number of photon in a particular mode is 
done by taking solution ^ and then calculating the mean 
value. This mean value should be computed with the 
initial density matrix, whose properties are discussed in 
the text. After straightforward calculations we get 



A: A 



(*)) 



167r 3 e ?iA 4 
x (f(q)\h(t)\ 2 



mm 2 ) 



(A5) 



where the function f(q) is defined in Eq. (fll]). This 
function is of the order of A^ 2 what may be easily demon- 
strated. Additionally the modulus of the function I/2I is 
much smaller than |/i|. Hence we are justified to neglect 
the contribution of function I/2I to our formula. In or- 
der to compute mean number of photons we have to sum 
over polarizations and integrate over frequencies. The 
first task may be performed with the help of well-known 

identity: X^a( 7 ^ ' ^k\) 2 = (•"■ — (^fc^d) 2 )- The informa- 
tion about the spectrum is hidden in the function fi(t). 
The integration over the frequencies may be easily per- 
formed if we assume the spectrum consist of very narrow 
lines. Since the pulse duration is long in comparison to 
the inverse of the optical frequencies: t 3> w£ ,u>q , we 
may approximate the spectrum to the form: 



I/1WI 5 



nA 2 t 



[5(uj k - lo ) + 5{u k - uj l )) 



(A6) 



As one can see the spectrum consists of two peaks cen- 
tered on the atomic and laser frequencies. The intensities 
of the two are the same, which is the result of the absence 
of the spontaneos emission during the short pulse. Now 
the result ( ^fj| ) is evident. The derivation of the Eq.jl^) 
describing the variance is done in the similar way as for 
the mean. 



APPENDIX B: NUMERICAL STRATEGIES 

The basic numerical task is to compute the functions 
f(q) and g(q) which enter into equations for the mean 
value (|l0|) and the variance ([l8]). For this purpose we 
express the mean value in the functions f(q) and g(q) by 
mean occupation numbers and higher moments of trap 
energy levels. We explain how it should be done for the 
function f(q), the derivation in the case of variance is 
similar. 

The summation in the function f(q) 



{9 H 9H'9 A 9m')Pfi,n' (q)Prn,,n> (Bl) 



n,n' ,m,m' 



is performed over four vector indices. We may exclude 
two of them because the mean of the four operators 
(9ft9ft' 9ih9m') is nonzero only if the operators couple into 
pairs. For four operators there are two possibilities of 
pairing, and we obtain 



f(q) = J2 (NnN^) (PnAfiPrnM-Q) + Ifc^P 



(B2) 



where Nft = g\\'gn- In derivation of ( p32| ) we use the iden- 
tity J2,n^MM)Pm,n'(-q) = 5n,n> that results directly 
from the completeness of the states. The Franck-Condon 
coefficients for the harmonic trap may be calculated ana- 
lytically. The result for three dimensional trap is simply 
the product of the ID Franck-Condon coefficients 



(n\e 



iqx I 



2^-2 1 a 



2 



(B3) 



where L™{x) denotes Laguerre polynomial and £ is the 
oscillator length. For simplicity we assume our trap 
to be spherically symmetric, therefore we may freely 
choose the orientation of the system of coordinates for 
the states in the trap. We choose the z-axis to be in 
the direction of the q vector. For such system of coordi- 
nates three dimensional Franck-Condon coefficients be- 
come {n\e l ^\rh) = (n z \e lqz \m z )S n:i;ym J nyymy . Since the 

mean values of the different combinations of Nft operators 
depends only on the energies of the states, and not on the 
distribution of the quantum numbers, we may split the 
summation into two steps: the sum over energies and the 
sum over degeneracies. The last one may be performed 
for all terms, and the result is a combination of Laguerre 
polynomials. For the microcanonical ensemble 



(A> 2 )-(A>o) (N E N E ,)(L E ( V )L E ,( V ) 



E<E' 



E>0 



+ E ((N 2 E )-(N E )-(N E N E )) S E ( V ) 



E>0 



N, (B4) 



where 77 = g 2 £ 2 /2 (the exponent in front of the square 
bracket is called the Debye factor) and the mean (N E N E ) 
is the correlation between different levels with the same 



energies. In the formula (B4) S^(x) denotes the sum 
including the Laguerre polynomials: 



(n + flt)! 



L a n (x) 2 (m+\-n). 



(B5) 



G 



The sum (B5) may be done analytically with use of 
the summation formula for Laguerre polynomials 8.974.1 
jl7j . The calculation leads to 

q«( T \ - ( m + 2 V- ( Ta +l T a ra+l ra ^_ ( m + X ) ! 



(m+a)! 



(m+a-1)) 



_J_ t a+2 t a — 1 jj_ t a+2 j a — 1 



1 



_ _ ra+i to 



1 



r a+1 r a 
-^m-l-^m-l 



(m+1)! 
(rn+a)! 



(B6) 



where we use shorthand notation = Lf n (x) and as- 
sume L^(i) = for m < 0. 
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